Renormalization factor and effective mass of the two-dimensional electron gas 
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We calculate the momentum distribution of the Fermi liquid phase of the homogeneous, two- 
dimensional electron gas. We show that, close to the Fermi surface, the momentum distribution 
of a finite system with A'^ electrons approaches its thermodynamic limit slowly, with leading order 
corrections scaling as N~^^^. These corrections dominate the extrapolation of the renormalization 
factor, Z, and the single particle elTective mass, m*, to the infinite system size. We show how 
convergence can be improved analytically. In the range 1 < Ts < 10, we get a lower renormalization 
factor Z and a higher elfective mass, m* > m, compared to the perturbative RPA values. 



PACS numbers: 

Landau's Fermi liquid theory [T] postulates a one-to- 
one mapping of low energy excitations of an interact- 
ing quantum system with that of an ideal Fermi gas via 
the distribution function of quasiparticles of momentum 
k. The resulting energy functional has been successfully 
applied to describe equilibrium and transport properties 
of quantum Fermi liquids, the most prominent are the 
electron gas and liquid ^He [2j [3]. However, quantita- 
tive microscopic calculations of its basic ingredients, the 
renormalization factor, Z, and the effective mass, to*, 
remain challenging. 

In this paper, we calculate these parameters for the 
2DEG (two-dimensional electron gas) using Quantum 
Monte Carlo (QMC) in the region 1 < Ts < 10, where 
= (7rna^)~^/^ is the Wigner-Seitz density parame- 
ter, n the density, and as — h^/lme^) the Bohr radius. 
Kwon et al.|3] made the first attempt to determine the 
Fermi liquid parameters of the two-dimensional electron 
gas using QMC, with results that differ from calculations 
based on other methods [5J[S1[7]. In particular, Kwon et 
al. found an effective mass smaller than the bare mass, 
e.g. TO* < TO at = 1. However, QMC calculations suf- 
fer from severe finite size effects since typical system sizes 
are limited to ~ 100 electrons. Here we show, that 
there is an extremely slow convergence of the effective 
mass and the renormalization factor to their thermody- 
namic limit values, with leading order corrections scaling 
as N~^^'^. A correct extrapolation to the infinite sized 
system leads to important qualitative and quantitative 
differences compared to previous calculations [H [7] which 
had assumed a (1/iV) extrapolation. We further use the 
knowledge of the analytical properties of the ground state 
wavefunction [S], to analytically estimate dominant and 
sub-dominant size effects which are important. 

Microscopically, the existence and characteristics of 
the Fermi surface of interacting fermions are directly re- 
lated to the renormalization factor Z at the Fermi sur- 



face |l9j. For a normal Fermi liquid, one expects a sharp 
discontinuity in the momentum distribution, nu, at the 
Fermi surface for each spin 

Ui, = Zi,0{kF - \k\) + gi, (!) 

where kp is the Fermi wavevector, 0{x) is a step-function, 
and (7k is a continuous function of momentum k. One ex- 
pects a linear excitation spectrum = h'^kp{k — kp)/m* 
close to the Fermi surface at |k| = kp, with to* the effec- 
tive mass. The goal is to determine the properties in the 
N —f oo limits based upon calculations on cells of N par- 
ticles with discrete values of and Z^ . A microscopic 
construction of the Landau energy functional is based 
on considering energy eigenstates which are adiabatically 
connected to the excited states of the non-interacting 
Fermi gas [7]. The effective mass is then explicitly given 
in terms of the energy difference, Ck = ^k^^^ ~ (^'o^ + m) i 
between single particle excitations of energy E^'^^ and 
momentum k, and the iV-particle ground state of energy 
Eq where /i is the chemical potential. 

Quantum Monte Carlo (QMC) methods provide the 
most accurate calculations of the ground state energy 
of the electron gas [101 E]. However, fermionic QMC 
calculations suffer from two major drawbacks, the fixed 
node approximation and finite size errors. For a normal 
Fermi liquid, the most precise results are obtained using a 
generalized Slater Jastrow form for the trial wavefunction 

m 

^-T oc i:i(R)e-^(^) (2) 

where R indicates a dependance on all particle coordi- 
nates. Antisymmetry is assured by a Slater determi- 
nant £'(R) = dctij e^^^'"^^^^^ of plane waves inside the 
Fermi sphere |kj| < kp using dressed quasiparticle coor- 
dinates qi(R) to account for many-body backflow effects, 
whereas the many-body Jastrow potential U (R) is sym- 
metric with respect to particle exchange and accounts for 
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FIG. 1: The momentum distribution of = 162 unpolarized 
electrons using grand-canonical twist-average VMC for densi- 
ties Ts = 1 — 10, analytically corrected for size effects around 
kp. The inset compares the uncorrected QMC data for dif- 
ferent system sizes between A'^ = 26 and A'^ = 162 with the 
size corrected distribution at = 10 ('W —> oo"). 



FIG. 2: The renormalization factor Z for — 10 estimated 
from the finite-size momentum distribution as a function of 
the inverse number of electrons, and the corresponding size 
corrected values. Dashed lines illustrate the size corrections 
of order N~^^'^ {^^^) for the uncorrected (corrected) data. 
The inset shows the corresponding values at ra = 1. 



the singularities in the interparticle potential at the coin- 
cidence points. Projector Monte Carlo methods (DMC) 
can be used to improve the wavefunction stochastically. 
Many ground-state properties have been successfully cal- 
culated using QMC, however, the situation is less clear 
concerning excited state properties [7]. 

The Slater determinant of the many-body wavefunc- 
tion, Eq. ([2]), directly connects the ground state of the 
interacting system with the non-interacting one: low- 
lying excitations are obtained by changing the "occupa- 
tion numbers" of the plane waves. The energy is there- 
fore a functional of the occupation numbers as postulated 
within Landau Fermi theory [3j. Whereas this energy 
functional certainly exists for any finite system, its exis- 
tence in the thermodynamic limit is non-trivial; a neces- 
sary condition is limTv^oo > 0; a central issue of this 
paper. 

We have performed Variational Quantum Monte Carlo 
calculations (VMC) of the 2DEG; the electrons inter- 
act with a 1/r potential and with a positive background 
charge. We have used a Slater- Jastrow backflow wave- 
functions (SJ-BF) with an analytical form for the both 
the Jastrow and backflow potentials [13]; all potentials 
are split in short and long-range contributions as de- 
scribed in [13]. For N = 58 electrons, the DMC ground 
state energies obtained are ^ 3mRy lower than previous 
calculations using numerically optimized forms [IF . Ex- 
cited states were formed by adding or subtracting orbitals 
in the determinant; the backflow and Jastrow forms [T3] 
are independent of the precise occupation of the Slater 
determinant. Since the trial function had no free pa- 
rameters, we can study size effects without re-optimizing 
parameters for different system sizes. 

First, we calculated the momentum distribution as ex- 
plained in Ref. [TT]. However, for systems in periodic 
boundary conditions, the momentum distribution is only 
given at discrete values k = 27r(nx -I- my) /L where n 



and m are integers and x, y are the unit vectors in the 
X and y direction, respectively. Using twisted boundary 
conditions with twist angles {9x^ + 0yy)2TT/ L for the trial 
wavcfunctions, we can obtain a momentum distribution 
for all values of k by varying the twist angle. In the limit 
of an infinite sized system, the Slater determinant of our 
trial wavefunction approaches a sharp Fermi surface, oc- 
cupying only wavevectors |k| < kp- For finite systems, 
the sharp behavior of the occupation numbers inside the 
Slater determinant is best described by working in the 
grand-canonical ensemble and for a given twist angle use 
only orbitals inside the Fermi sphere. This leads to a 
varying particle number as a function of the twist angle. 
As described in Ref. [8J, the translational invariance of 
the ground state wavefunction allows us to define pock- 
ets inside of which the wavefunction transforms trivially 
- any change of the twist angle inside a pocket reduces to 
a change of the total center of mass moment accounted 
for by a simple phase factor; only a single QMC calcula- 
tion is needed for each pocket. As shown in Fig. [T| the 
renormalization factor quantifying the jump in the mo- 
mentum distribution at kp can be read-off precisely for 
any finite system. However, strong size effects around 
the Fermi surface are still evident. 

We can analyze size-effects directly using the analytical 
form of the SJ-BF trial wavefunction. The momentum 
distribution is obtained by displacing one particle Vj a 
distance r: 

< - ,3, 

where (...) denotes averaging over |^'t(R)P and over a 
uniform distribution for r. The change of the Jastrow 
factor in Eq. ^ writes 
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where pq = J^j e^^'^^ ■ As described in Ref. |S] , the most 
important finite size-effects can be understood as an inte- 
gration error by analytical continuation of the finite-size 
(periodic) wavefunction to an infinite system where the 
estimator in Eq. (|3| would contain the following change 
in the Jastrow factor 



(27r)2 



u,[e^-^'p^^-l\ [e:^--l]. (5) 



The finite-size correction is then dominated by the non- 
analyticity of the integrand at (7 = 
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FIG. 3: The inverse effective mass ra/m* for rs = 10 as a func- 
[e*'* '^^ (0_q — 1] [e*'''''" — 1] , tion of A'^"^, together with the corresponding size corrected 

values. Dashed lines illustrate the expected size corrections. 



and we can calculate the leading order size corrections, 
(5nk = — , by expanding , Eq. ([s]) , up to second 
order in SUoo — SUn- Neglecting mode-coupling terms, 
we get 



k+q ~ "-k . 



T/i 

yr/L ,2 



^UqUqi [1 - S{q) - S{q)\ 



+q "k+q'. 



where 



5{q) = K (1 - S{q)) 



nu1S{q)\ 



(6) 



(7) 



Equation ([6]) expresses size corrections in terms of the 
long-wavelength limits of the Jastrow potential and the 
structure factor, S{q). For the homogenous electron gas, 
in the limit q —t Q, we have: 



2nuq ~ -1 + [1 + (2nWq/eg)]i/2 
S{q) ~ [2nu, + l/So{q)]-^ 



(8) 



where Vq — 2TTe^/q, £q — h^q^/2ra, and So{q) is the 
structure factor of the non-interacting Fermi gas. 

As the momentum distribution of a Fermi liquid, 
Eq. ([1]), is smooth everywhere away from the Fermi sur- 
face, leading order corrections are restricted to a small 
region around kp, where we can write 



In Figure [1] we show the size-corrected momentum distri- 
bution for different densities between rs = 1 and rs = 10 
using Eq. (|9|. Close to kp, size effects lead to important 
qualitative and quantitative changes. 

The renormalization factor, Z, can be read-off di- 
rectly from the jump of the momentum distribution at 



the Fermi surface, Zj^ 



n^„ , — ?T-^_[_£, and its size 



corrected value may therefore be read of directly from 



Fig. [T] For a precise evaluation of Z in the thermo- 
dynamic limit, we have studied the extrapolation sep- 
arately. From Eq. ([6| , one can show that size-corrections 
of Z can be written as 



An = 



dh 



-.IL (2^)- 



(10) 



which includes the main sub-leading order corrections. 
Using the analytical forms, Eq. Q and Eq. ([s]), the lead- 
ing order corrections are 



2 \ -1/4 
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(11) 



The asymptotic form, Eq. (10 1 with Eq. (11 1, shows that 
actual QMC calculations with typically N ~ 10^ elec- 
trons suffer from very strong size effects. Obscured by 
the intrinsic noise of QMC calculations, pure numerical 
analysis of the data might suggest convergence to values 
far off the exact value in the thermodynamic limit. 

In Figure [2] we compare the bare data for = 1 
and rs = 10 with their size corrected values. Whereas 
the bare data are in reasonable agreement with previ- 
ous QMC results [3 [TT], a numerical extrapolation of 
the uncorrected data strongly depends on assumptions 
on the asymptotic scaling form, as size corrections over- 
whelmingly dominate the calculation oi Z. In order to 



go beyond leading order, we have directly used Eq. (10 1 
together with the asymptotic forms, Eq. ([?]) and Eq. |8f^, 
to correct our bare data analytically. As can be seen from 
the figure, the size corrected values drastically reduce size 
effects, as expected. More important, in contrast to the 
uncorrected data, the extrapolation of the size corrected 
values is not sensitive to assumptions on the remaining 
corrections for densities rs > 3. Approaching the high 
density region rs ^ 1, the thermodynamic limit extrapo- 
lation is getting more difficult, since the asymptotic ex- 
pansion is singular in the limit Ts 0. In table|T]we have 
summarized our results for the renormalization factor. 
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1 


3 


5 


10 


-EvMC 


-0.4179(2) 


-0.4223(1) 


-0.2975(1) 


-0.16952(1) 


-Edmc 


-0.4206(2) 




-0.2991(1) 


-0.17070(1) 


Zvmc 


0.62(4) 


0.34(3) 


0.22(2) 


0.090(4) 


m*/mvMC 


1.26(7) 


1.39(8) 


1.54(7) 


1.72(9) 


^RPA 


0.66 


0.44 


0.34 


0.24 


m*/mRPA 


1.02 


1.12 


1.16 


1.21 



TABLE I: Energies per particle (in Ry), -Bvmc, and -Edmc, 
the renormalization factor Z and the effective mass m* /m 
extrapolated to the thermodynamic limit (both within VMC, 
and from perturbative RPA calculations [5]). Values in () are 
standard errors in the last decimal place. 



Size corrections of the momentum distribution induces 
size corrections for the total kinetic energy which can 
be shown to coincide with the two-dimensional analog of 
Ref. [5]. In two dimensions, leading order size corrections 
of the kinetic and potential energy per particle scales as 
7V~^/4 in the Fermi hquid phase. We have added VMC 
and DMC energies of the size-extrapolated values of the 
energy per particle in table |T]. 

Since our class of wave function have Z > Q, the sin- 
gle particle excitation spectrum should be dominated by 
quasiparticle excitations with an effective mass m* . We 
have calculated the effective mass by adding an electron 
with momentum p with |p| > kp to the ground state. 
The effective mass of an excited state has been deter- 
mined assuming an expansion of the self energy in powers 
of p— kp, leading to ImSplfi' = — kp + 2kp(m/m* — 
l)(p — kp) in the vicinity of kp. 

Again, the proper treatment of size effects is essential 
to extrapolate to the thermodynamic limit. The addi- 
tional electron at momentum p will induce size correc- 
tions in the momentum distribution which can be esti- 
mated as before. The resulting additional finite size error 



in the total kinetic energy, 6Tp , due to the excitation of 
momentum p, is then given by 



P 2to P I- 



1] 



(12) 



We see that size-corrections of the effective mass are in- 
trinsically related to those of the renormalization factor, 
Z, leading to a similar asymptotic scaling law, iV^^/^. 
Potential energy corrections are independent of p in lead- 



ing order, and Eq. ( 12 1 dominates finite size corrections 
for m/m* . Note, that the renormalization factor can also 
be obtained from analyzing the finite-size error of effec- 
tive mass calculations without explicit calculations of the 
momentum distribution. 

From figure |3] we see that size effects play a similar 
important role for determining m* as they do for deter- 
mining Z. In particular, for high densities, size effects 
qualitatively change the conclusion of previous calcula- 
tions H]: whereas, in agreement with f?^ all bare data 



indicate an effective mass smaller than the bare mass for 
N ^ 100, in the thermodynamic limit the effective mass 
is increased, as predicted by perturbative RPA calcula- 
tions [F,^. 

Calculations based on many-body perturbation theory 
going beyond the perturbative RPA approximation have 
been suggested. However, based on different approxima- 
tions, these predictions may lead to an enhancement or 
depression of Z (or the effective mass) [51 1^ and it is dif- 
ficult to estimate reliable the validity of the underlying 
approximations. Our VMC results for Z are always below 
the corresponding values of the perturbative RPA calcu- 
lations, whereas we predict a higher effective mass m* /m. 
Our calculations therefore support improved RPA calcu- 
lations based on many-body local field theory including 
charge- and spin-density fluctuations as proposed in [6 . 

This research was supported by NSF DMR04-04853, 
IDRIS Computers, and the ACI "Desordre et Interac- 
tions Coulombicnncs" , and facilitated by the Project de 
Collaboration CNRS/UIUC. M.H. thanks S. Moroni, C. 
Pierleoni, R. Chitra and A. Pasturel for discussions. 



[1 
[2 

[3; 
[4; 

[5 
[6 

[7: 

[8 
[9 

[10 

[11 

[12 
[13 
[14 

[15 



L. D. Landau, JETP. (USSR) 30, 1058 (1956). 
D. Pines and P. Nozieres, The Theory of Quantum Liq- 
uids, Addison- Wesley, 1989. 

G. Baym and C. Pethick Landau Fermi Liquid Theory, 
Wiley- VCH, 2004. 

Y. Kwon, D.M. Ceperley, and R. M. Martin, Phys. Rev. 
B 50, 1684 (1994); Phys. Rev. B 53, 7376 (1996). 

H. -J. Schulze, P. Schuck, and N. Van Giai, Phys. Rev. B 
61, 8026 (2000). 

R. Asgari, B. Davoudi, M. Polini, G. F. Giuhani, M.P. 
Tosi, and G. Vignale, Phys. Rev. B 71, 045323 (2005). 
G. F. Giuliani and G. Vignale, Quantum Theory of the 
Electron Ltqmd, Cambridge University Press, Cambridge 
(2005). 

S. Chiesa, D.M. Ceperley, R.M. Martin, and M. Holz- 
mann, Phys. Rev. Lett. 97, 076404 (2006). 
V.M. Galursky and A.B. Migdal JETP (USSR) 34, 139 
(1958) - Soviet Phys. JETP 7, 96 (1958); J.M. Luttinger, 
Phys. Rev. 119, 1153 (1960). 

D.M. Ceperley and B. Alder, Phys. Rev. Lett. 45, 566 
(1980). 

B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, 5005 
(1989). 

M. Holzmann, B. Bernu, and D.M. Ceperley, Phys. Rev. 
B 74, 104510 (2006). 

M. Holzmann, D.M. Ceperley, C. Pierleoni, and K. Esler, 
Phys. Rev. E 68, 046707 (2003). 

V. Natoh and D.M. Ceperley, J. Comput. Phys. 117, 171 
(1995); M. Holzmann and B. Bernu, J. Comput. Phys. 
206, 111 (2005). 

Y. Kwon, D.M. Ceperley, and R. M. Martin, Phys. Rev. 
B 48, 12037 (1993). 



